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ABSTRACT 



In the double neutron star merger or neutron star-black hole merger model 
for short GRBs, the outflow launched might be mildly magnetized and neutron 
rich. The magnetized neutron-rich outflow will be accelerated by the magnetic 
and thermal pressure and may form a two component jet finally, as suggested by 
Vlahakis, Peng & Konigl (2003). We show in this work that such a two component 
jet model could well reproduce the multi-wavelength afterglow lightcurves, in 
particular the X-ray flat segment, of short GRB 051221A. In this model, the 
central engine need not to be active much longer than the prompt 7— ray emission. 

Subject headings: gamma-rays: bursts - ISM: jets and outflows - binaries: close 
- stars: neutron 



1. Introduction 



Kouveliotou et al. (1993) showed that the Gamma-ray Bursts (GRBs) can be divided 
into two distinguished groups: long GRBs (LGRBs) with duration > 2 sec and short GRBs 
(SGRBs) with duration < 2 sec. Fruitful afterglow data have been collected for LGRBs 
and support the collapsar model (see Woosley & Bloom 2006 for a recent review). The 
afterg low data of SGRBs are relatively rare and have not been reliably detected before May 
2005 (IGehrels et al. 20051 ). The l eading model for SGRBs is the merger of neutron stars or 
neutron star-black hole binaries (jEichler et al. 19891 : iNarayan et al. 1992l ). One attractive 
feature of such a compact star binaries merger model is that the accretion disk formed is 
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too small to give rise to a long burst, provided that the central remnant has collapsed to a 
black hole immediately (e.g., Narayan et al. 2001: Rosswog et al. 2002; Aloy et al. 2005; cf. 
Kluzniak & Ruderman 1998; Rosswog et al. 2003). Recently, this view of SGRB origin is 
supported by observations (see Meszaros 2006 for a recent review) : a) For some short bursts, 
the host is early type galax y, with a stellar population older than ~ 1 Gyr. The overall star 
formation rate is very low (IFox et al. 20051 ; iBarthelmy et al. 20051 ). b) Deep search of an 
underlying supernova has resulted in null results (e.g. Hjorth et al. 2005). c) The typical 
number density of the medium surrounding the SGRB progenitors is low to 10 -3 cm -3 and 
a stellar wind profile is disfavored. 



While the compact objects binary merger model for SGRBs is largely consistent with 
the ongoing afterglow observations, some lightcurve data, for example, the optica l or X-ray 
flares in GRB050709 (Fox et al. 2005) and G RB050724 faarthelmv et al. 2005h. and the 
X-ray afterglow flat segment in GRB051221A (jSoderberg et al. 20061 ; iBurrows et al. 20061 ). 
have not been well interpreted. One sp eculation is that the central engine still works after 



the cease of the prompt 7— ray em ission (IFan et al. 20051 ; IZhang et al. 20061 ; iDai et al. 2006 



Gao fc Fan 20061 ; IFan &: Xu 20061 ) and may involve some magnetic processes. For example, 
Fan & Xu (2006) suggested a magnetar wind energy injection model for the X-ray afterglow 
flat segment lasting ~ 10 4 seconds in GRB 051221A. However, it is still unclear that whether 
the supermassive or hypermassive magnetar formed in the double neutron s tar coalescence 
(Price & Rosswog 2006) could survive for such a long time (IDuez et al. 20061 ). In this Letter 
we suggest a two component jet model to fit the multi-wavelength afterglow data of GRB 
051221A. In this model, the central engine need not to be active much longer than the 
prompt 7— ray emission. We discuss the possible physical scenario giving rise to such a 
two component jet in section 2. We fit to the afterglow data numerically in section 3. We 
summarize our result with some discussions in section 4. 



2. A two component jet model 



In the double neutron stars or neutron star-black hole merger scenario, the ou tflow is 
very likely to be neutron rich and might also be magnetized (IPrice fc Rosswog 20061 ). For a 
magnetized neutron-rich outflow, the accel eration process is significantly different from that 
of the standard fireball (IPiran et al. 1993T ). As shown in Vlahakis, Peng & Konigl (2003), 
the neutron-to-proton ratio in disk-fed outflows might be as high as ~ 30 and the neutrons 
can decouple at a Lorentz factor T ~ tens, while the protons continue to be accelerated to 
a T ~ hundreds and is collimated to a narro wer angle by the electromagnetic force. As a 
result, a two component jet might be formed (I Vlahakis et al. 20031 ; iPeng et al. 20051 ). We 
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use such a structured jet profile to fit the afterglow of short GRB051221A in next section. 



3. Fit to the afterglow light curve of GRB 051221A 



3.1. The observation 



GRB051221A located by Swift is one of the few S GRBs whose afterglows were well 
detected in multi-bands. It s duration is T 3n = 1 .4 ± 0.2s (jCummings et al. 20051 ). the peak 
energy is E p = A02^jlkeV (IBurrows et al. 2005), th e gamma-ray (20keV-2MeV) fluence is 



about 3.2 x 10 6 erg ■ cm 2 (jGolenetskii et al. 20051 ). With the optical spectrometric mea- 



sured redshift z = 0.5464 JSoderberg et al. 2006h . the total isotropic gamma-ray (20keV- 
2MeV) energy of GRB051221A is about 2.4 x 10 51 erg, make it th e most fluent SGRB and 
about 6 — 35 times higher than the others (jSoderberg et al. 20061 ). The Swift XRT began 
observing GRB051221A 88 seco nds after the burst and found the bright, rapidly fading X-ray 
afterglow (IBurrows et al. 20051 ). The X-ray afterglow early decay index is «xi = 1-I6l j7 
(100-2200s) and late a X3 = 1.09±g;Jg (3 x 10 4 - 3.6 x 10 5 s), between them is a slow de- 
cay decline «X2 = 0.041q;|[, an d the average photon spectral index of these three segments 
is T x = 2.1 ± 0.2(Burrows et al. 2006). The Chandra observations show a steep decay 
a X 4 = 

is r x 



1.93+ q'iq a fter ~ 3.6 x 10 5 s, an d the average photon spectral index of this segment 



-0.19 
-0.29 



1.9415$ ( IBurrows et al. 20061 ). The GMOS start optic al observation in r' ban d 2.8 



hours after burst and get an optical decline a opt = 0.92 ± 0.04 (jSoderberg et al. 20061 ). 



The spectral analysis of the optical-to- X-ray afterglow data at ~ 1 day after the trigge r 
of the burst suggests a cooling frequency v c ~ 2 ± 1 x 10 17 Hz (jSoderberg et al. 20061 ). 
With the early X-ray decay index (100-2200s), we have an electron index p = 1 — 4oxi/3 = 
2.55+o;23, note that v c oc t~ 1/2 is over the XRT band. And later (3 x 10 4 - 3.6 x 10 5 s), 
when v c is crossing the X-ray band, the decay index smoothly changes from (3 — 3p)/4 
to (2 — 3p)/4. The average decay index «X3 is between these two values, that we have 



2.12 



hO.13 
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(2 - 4a X3 )/3 < p < 1 - 4a X3 /3 = 2.45; 



hO.13 
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Or by the optical decay index at the 



same time, we have p — 1 — 4a opt /3 = 2.23 ± 0.05 (u c is still above the optical band now). 
It is interesting to note that the early and later electron distribution indexes p seem to be 
much different (see also Burrows et al. 2006 for a similar argument). In previous numerical 
fit to the afterglow data of GRB 05 1221 A, a constant p can not well reproduce both the 
early and the late X-ray declines, as shown in Figure 3 of Burrows et al. (2006) and Figure 
1 of Fan & Xu (2006). In our two component jet model, the shock parameters (e e , €b, p) 
might be different for the narrow and the wide components. The early afterglow emission is 
mainly contributed by the narrow jet component while the late afterglow is dominated by 
the forward shock emission of the wide jet component. So the values of p for the early and 
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later afterglow data could be different. 



3.2. Numerical method 



We interpret the apparent flattening in the X-ray light curve being caused by the emer- 
gence of the forward shock emission of the wide jet component. The code here to fit the 
multi-band lightcurves has been used in Yan et al. (2007), in which the outflow dynamical 
evolution is governed by a set of differential equations. 

Reverse shock phase. In this case, the energy conservation of the whole system that 
contains four regions (1) the unshocked medium, (2) the shocked medium, (3) the shocked 
ejecta, and (4) the unshocked ejecta (reads Huang et al. 2000; Yan et al. 2007) 



d{E 2 + E 3 + E A ) = -e 2l2 dU 2 - £ 3 73^3, (1) 

where E- Y is kinetic energy measured in observer frame and they are governed by E 2 = 
(72 -l)m 2 c 2 + (1- £2)72^2, £3 = (73 -l)m 3 c 2 + (1-63)73^3 and E A = (77 - l)(m cj - m 3 )c 2 . 
Where U 2 = (72 — l)m 2 c 2 and U 3 = (734 — l)m 3 c 2 are internal energy measured in the 
comoving frame. Note that in this subsection, the subscripts i = 1 — 4 represent the i— th 
region, respectively. E\ is the radiation fraction of internal energy, 7; are the velocities (in 
units of c) and the corresponding Lorentz factors measured in observer frame, respectively; 
rrii is the rest mass of the material, 734 is the Lorentz factor of region 3 relative to region 4. 

The mass of the medium swept by the forward shock is 

dm 2 = A7rR 2 nim p dR, (2) 

where n\ is the proton number density measured in the comoving frame, m p is the proton's 
rest mass and R is the radius of forward shock front to the central engine. 

The mass of the medium swept by the reverse shock is 

dm 3 = 4nR 2 r}n4m p (P4 — P^s)dR, (3) 

where 77 is the initial Lorentz factor of outflow, and /3rs ~ (73^3/^3 — 74n 4 /?4)/(73n 3 — 7477-4) 
( jSari fc Piran 19951 ; iFan et al. 20041 ) is the velocity of reverse shock. 

With eq.(l-3) and assuming that 72 = 73, we have the following equation to describe 
the dynamical evolution of forward shock 

j A D 2 (72 - l)nimp + (72734 - V)v n <tm p (fl4 ~ Ais) m 
«T2 = -4:TtR , (4) 
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where / = m2+m 3 + (l-e2)(2j2-l)m2 + (l-e 3 )(j 3 4-l)m 3 + (l-e 3 )-f2m 3 [r](l-p 2 P4) 



After the cease of the reverse shock. In t his case, the dynam ical evolution of the 
forward shock has been discussed in great detail (IHuang et al. 20001 ) . The reverse shock 
electrons cool adiabatically. We take the analytical approach of Sari & Piran (1999) to 
calculate the emission of those electrons. 

In the afterglow lightcurve calculation (see Yan et al. 2007 for details), the cooling of 
the shoc ked electrons due to both synchrotron and inverse Compton has been taken into 
account (jSari et al. 19981 ; IWei &: Lu 19981 ; ISari fc Esin 200 ll ). The synchrotron radiation of 
the forward shock electrons on the "equal arrivin g surface" (on whic h the emission reaches 
us at the same time) has been integrated strictly (IHuang et al. 20001 ) . 



3.3. Fit to the afterglow data of GRB051221A 

Using the code described above, we fit the afterglow of GRB051221a. Here we consider 
two c ases. First, we assume that the shock physical parameter s of two components are the 



same flBerger et al. 20031 : ISheth et al. 2003 



Huang et al. 2004 ). The parameters taken into 



account are as follow: for the narrow component, the isotropic energy -Ek.iso = 3.2 x 10 51 
ergs, the initial Lorentz factor 77 = 500, and the half opening angle Bj = 0.03; for the 
wide component, -Ek,iso = 9.5 x 10 51 ergs, r\ = 50, and Bj = 0.1. For both components, 
e e = 0.3, €b = 0.01, and p = 2.2. However, from Fig.l we note that in this case the fit 
is not very well, the reason is that we take the same value of p for the two components 
(similar conclusion can be found in Burrows et al. 2006 and Fan & Xu 2006). So we take 
a different value of p = 2.8 for the narrow component and find that the numerical results 
can well reproduces the afterglow of GRB051221A. In this case the ratio of total energy of 
narrow and wide component is 1 : 30, which is very different from 3 : 7, the ratio found in 
Vlahakis et al.'s simulation (2003). Second, we assume that for the narrow and the wide 
components, the shock parameters are different. This assumption is somewhat unusual. For 
the pr e- Swift GRBs, the best fitt e d micro physical parameters are different from burst to 



burst (IPanaitescu fc Kumar 200ll : lYost et al. 20031 ). For the Swift GRBs, the situation is 



less clear because the afterglows are usually too peculiar to be reasonably reproduced in 
the standard fireball model (e.g., Fan & Piran 2006). Figure 2 Shows our numerical fit to 
the X-ray and optical afterglows with the following parameters: for the narrow component, 
the isotropic energy -Ek,iso = 3.2 x 10 52 ergs, the initial Lorentz factor rj = 500, Bj = 0.03, 
e e = 0.09, €b = 0.003, and p = 2.8; for the wide component, -Ek,iso = 9.5 x 10 51 ergs, r] = 50, 
Bj = 0.1, e e = 0.3, e B = 0.01, and p = 2.2. 



- 6- 



4. Summary and Discussion 



In this work we show that the multi-band afterglows of GRB051221A could be well 
fitted by a two component jet model. Such a t wo component jet m ight be formed if the 
disc- fed outflow is neutron rich and magnetized (IVlahakis et al. 20031 ). Our good fit to the 
afterglow data in turn suggests that the initial ejecta launched in the double neutron stars 
merger or neutron star-black hole merger may be neutron rich and magnetized. However, 
we point out that any model consisting of a fast and a slower jet can also provide a good fit 
to the observed data for 051221A. This and some other possibilities are described below. 

There are several possible ways to produce a slower jet: 1) Neutrino annihilation and 
magnetic field were considered as the ene rgy source for SGRB in NS binary merger scenario 
(jEichler et al. 19891 ; iNarayan et al. 19921 ). Maybe they both operate, but only one mecha- 
nism produced the GRB prompt emission wh i le the other contribute only to late afterglow 
as a second component (IRosswog et al. 20031 ; iPrice fc Rosswog 20061). 2) A thick d isc acts 
as high-density wall was considered to collimate GRB outflow ( IRosswog et al. 20031 ). It was 
thought if mixing instabilities at the walls spoil the jet with baryons, it would result in a 
non-relativistic "dirty" fireball (IRosswog et al. 20031 ). We consider that if the disc changes 
to unstable when the GRB outflow launch, the consequent outflow might be more baryon 
loaded and less collimated, maybe as a second component. It should be noted that the 
validity of these processes are quite uncertain and more detailed investigation are needed 
before they can be used for afterglow calculations. 

We would like to point out that the two component jet model is not unique to account 
for the afterglow data of GRB051221A. One possibility is that the GRB outflow might have 
a range of bulk Lorentz fa ctors and the inner part s moving with a bulk Lorentz factor ~ tens 
carry most of the energy (IRees fc Meszaros 19981 ) . As the outer faster parts got decelerated 
by the external medium, the inner parts would catch up with the decelerating forward shock 
and increase its kinetic energy. For the particular GRB 05 1221 A, an energy-lorentz factor 
distribution E(> V ~ 20) tx r~ 4 5 is needed to reproduce the X-ray afterglow flat segment 



(jSoderberg et al. 20061 ; iBurrows et al. 20061 ). But the physical process pumping such kind 
of energy injection, in particular in the SGRB scenario, is not clear yet. The other possibility 
is that the central engine is active much longer than the prompt 7— ray emission phase. For 
example, Fan & Xu (2006) proposed a magnetar wind energy injection model to account for 
the long term X-ray flat segment. However, it is still unclear that whether the supermassive 
or hypermassive magnetar formed in the double neutron star coalescence could survive for 
such a long time. In principle, there could be a method to distinguish the two component 
jet model from these two possibilities — In the two component jet model, the first jet break 
might be observed when the edge of the narrow component enters our line of sight. The lack 
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of a detection of the first jet break renders the energy injection model possibly. The future 
fruitful muti-wavelength afterglow data may help us to pin down the underlying physical 
processes. 
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Fig. 1. — Fit to the multi-band afterglow of GRB051221A. Circles and squares are optical(r' 
band) and X-ray (0.3-10keV) observations. The solid lines are our calculation. To make a 
comparison, we also take a different p = 2.8 for the narrow jet, as the dashed lines shown. 
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Fig. 2. — Fit to the multi-band afterglow of GRB051221A. Circles and squares are optical(r' 
band) and X-ray (0.3-10keV) observations. The dashed, dotted and solid lines are narrow, 
wide jet and combined calculation. 



